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Abstract: For some time, the Draco dwarf spheroidal galaxy has garnered interest as a pos- 
sible source for the indirect detection of dark matter. Its large mass-to-light ratio and relative 
proximity to the Earth provide favorable conditions for the production of detectable gamma 
rays from dark matter self-annihilation in its core. The Solar Tower Atmospheric Cherenkov 
Effect Experiment (STACEE) is an air-shower Cherenkov telescope located in Albuquerque, 
NM capable of detecting gamma rays at energies above 100 GeV. We present the results of 
the STACEE observations of Draco during the 2005-2006 observing season totaling 10 hours 
of livetime after cuts. 



Introduction 

Dark matter is thought to be an important 
component of the Universe and research into 
its nature is actively pursued using a variety of 
techniques. Dark matter may be weakly inter- 
acting massive particles (WIMPs) which would 
tend to accumulate at the bottom of gravita- 
tional potential wells, such as galaxies, where 
they could undergo self-annihilation processes. 
Depending on WIMP mass and branching ra- 
tios, a measurable flux of high energy gamma 
rays could result. 

The Draco dwarf spheroidal galaxy has long 
garnered interest as a potential source of con- 
centrated dark matter [1]. Draco has one of 
the highest known mass-to- light ratios (M/L), 



perhaps as high as SOOMq /Lq [2] . Current ob- 
servations are consistent with a cuspy density 
profile [3], which would enhance the gamma- 
ray production rate. Furthermore, since Draco 
is a satellite of the Milky Way, its relative prox- 
imity to the Earth {d ^ 75 kpc) [4] might allow 
for the detection of such gamma rays. 

STACEE Observations of Draco 

The Solar Tower Atmospheric Cherenkov Ef- 
fect Experiment (STACEE) is a gamma- 
ray telescope operating at the National So- 
lar Thermal Test Facility (NSTTF) in Al- 
buquerque, NM. STACEE is a wavefront- 
sampling Cherenkov telescope which uses 64 
of the mirrors in the NSTTF heliostat array 
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ON events 


OFF events 


Excess 


Significance 


After Time Cuts 


177498 


177273 


225 


+0.39(7 


+ grid ratio Cut 


3094 


3120 


-26 


-0.33(7 



Table 1: Data summary of STACEE observations of Draco during the 2005-2006 observing season, 
representing approximately 3.67 x 10* s of livetinie. 




Figure 1: Effective area curves for STACEE 
observations of Draco. The blue (solid) line 
represents the STACEE effective area with- 
out cuts, the red (dashed) hne represents the 
STACEE effective area including the grid-ratio 
cut. 



for a total of ~ 2400 of collecting sur- 
face. Cherenkov light from gamma-induced air 
showers is reflected off the heliostats onto sec- 
ondary mirrors on a tower on the south side 
of the field. These secondaries focus the light 
onto photomultiplier tubes (PMTs) in such a 
way that each PMT sees the light from a single 
heliostat. Pulses from the PMTs are split, with 
one copy discriminated and used in the forma- 
tion of a trigger and the other digitized us- 
ing a 1 GS/s FADC. The trigger selects show- 
ers that deposit light evenly over the heliostat 
field, which favors those showers initiated by 
gamma rays over those resulting from charged 
cosmic rays, the most important background 
for the STACEE experiment. For a more com- 
plete description of the STACEE experiment, 
see [5]. 

The basic unit of observation for STACEE is 
the "ON-OFF" pair; 28 minutes on-source and 



28 minutes off-source. Both observations view 
the same path across the sky in local coordi- 
nates (altitude and azimuth), but separated by 
30 minutes in celestial coordinates (right as- 
cension). The off-source observation allows for 
a measurement of the local background condi- 
tions. We measure the significance of a mea- 
surement as 

ON-OFF S 

a = , : ^ , (1) 

VOWTOFF V2B 

where S is the signal and B is the background. 

STACEE observations of Draco total 35 "ON- 
OFF" pairs, of which approximately 10 hours 
of livetime remain after excluding periods with 
bad weather and known technical difficulties. 
Our data set is summarized in Table 1. 

Data Analysis 

Our raw background trigger rate from cosmic 
rays is approximately 5 Hz. In order to reduce 
this, we perform a grid-ratio cut which prefer- 
entially removes hadron-induced showers. This 
technique has been used successfully by the 
CELESTE experiment [6] and our implemen- 
tation is described in more detail in [7]. A 
simplistic description of the technique is that 
the "smoothness" of a shower is measured by 
the height-to-width ratio (H/W) of the sum 
of pulses from all 64 channels in the detector. 
This quantity depends on the relative timing 
of each FADC trace, which depends on the as- 
sumed impact point of the shower core (i.e., the 
extrapolated shower axis). The grid-ratio cut 
is based on how sharply peaked the H/W dis- 
tribution as a function of assumed core position 
is. Gamma-ray showers, smooth and symmet- 
ric, are expected to produce narrower H/W 
distributions than hadronic showers, which re- 
sult in broader, clumpier deposits of Cherenkov 
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light. Applied to our 2002-2004 Crab data, 
the grid-ratio cut improves the detection sig- 
nificance from 4.8cr to 8.1cr[8]. 

As seen in Tabic 1, wc do not detect an ex- 
cess gamma-ray signal from Draco in our data 
set. We derive an upper limit for the flux from 
Draco given a measure of our detector response 
to a candidate source spectrum. We discuss 
two possible source spectra, an E~^ '^ power 
law (suggested by the gamma-ray flux from the 
galactic center [9]) and a candidate dark matter 
spectrum taken from Tyler [1], with an energy- 
dependent shape given by 

'^("<^^= . (.^-5 +0.00014) ' 

where x = E/m-^c?. This gives a sharp cutofi^ 
at the energy corresponding to the candidate 
WIMP mass. 

Power Law Spectrum: Figure 1 shows ef- 
fective area curves for STACEE observations of 

Draco. Gamma-ray showers were simulated us- 
ing the Monte Carlo chain of the CORSIKA air 
shower simulation together with our own opti- 
cal ray-tracing model for the heliostats, sec- 
ondaries, and PMTs, and a simulation of the 
electronics [8, 10]. The effective area is given 
by the product of the probability that a shower 
triggers our detector and the area over which 
the simulated showers were generated. 

STACEE has an energy-dependent response 
which means the sensitivity to a given source 
depends on its energy spectrum. Figure 2 
shows the result of convolving the effective area 
curves with candidate spectra. As is custom- 
ary, we define the energy thresholds of our mea- 
surements as the peak of these curves. 

The flux limit is defined by 



NuL= T 



A{E)^{E)dE (3) 



where Nul is given by the 95% upper limit 
of the excess Nqn — Nqff, T is the livetime, 
and A{E) is the effective area. The differential 
fiux, ^{E) — C(j){E), contains a normalization 
constant with units of [cm~^ GeV~^]. 
For the data given in Table 1, including the 
grid-ratio cut, Nul = 138 and the result- 
ing upper limit is < 4 x 10~^ ( cfi/) 



7 cm~'^ GeV~^ at a characteristic energy 
of 220 GeV. Figure 3 shows a comparison of 
this limit with the published upper limit of the 
Whipple collaboration[ll] and our own Crab 
spectrum. 

Tyler Spectrum: Tyler provides a prescrip- 
tion for converting a flux limit into a cross- 
section limit (< av >xx) assuming a spherical 
isothermal halo model where the mass density 
is given by phaio = Ar~^ . We avoid a diver- 
gence at the center by including a constant- 
density core physically motivated by an equi- 
librium between infalling particles and annihi- 
lation; 



1/2 

= Rext < crv >xx 



Phali 



1/4 



(4) 



where Rext is the outer radius of Draco and 
< (TV >xx is the expectation value of the self- 
annihilation rate, given by the product of the 
cross-section and the velocity of the WIMPs in 
the halo. We then substitute this into our flux 
equation 



4^2 



< av >xx ME) (5) 



Sd'^Rr, 

and solve for the self-annihilation rate. The 
resulting limits are shown in Figure 4. 

Conclusions 

STACEE does not detect a significant sig- 
nal from Draco, and sets upper limits on 
cross-sections for WIMP with rest-mass energy 
greater than about 150 GeV. 
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Figure 2: Response Functions for STACEE observations of Draco given two candidate spectra. 
The left figure corresponds to an spectrum, and the right figure corresponds to a Tyler 

spectrum (Eq. 2) for an example 300 GeVjc? WIMP. The blue (sohd) curves represent the STACEE 
sensitivity without cuts, the red (dashed) curves include the grid-ratio cut. 




Figure 3: STACEE Flux limits for a £;~2.2 
spectrum as applied to Draco (blue) . For com- 
parison, also shown is the energy spectrum 
of the Crab Nebula (green) as measured by 
STACEE during 2003-2005 which is weU fit by 
the form: 4i = 1-2 x lO"'^ £;"2.23 
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